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ABSTRACT 

The problem of quantizing a class of two-dimensional integrable quantum field theories 
is considered. The classical equations of the theory are the complex sl(n) affine Toda 
equations which admit soliton solutions with real masses. The classical scattering theory of 
the solitons is developed using Hirota's solution techniques. A form for the soliton S-matrix 
is proposed based on the constraints of S-matrix theory, integrability and the requirement 
that the semi-classical limit is consistent with the semi-classical WKB quantization of the 
classical scattering theory. The proposed S'-matrix is an intertwiner of the quantum group 
associated to sl(n), where the deformation parameter is a function of the coupling constant. 
It is further shown that the S'-matrix describes a non-unitary theory, which reflects the fact 
that the classical Hamiltonian is complex. The spectrum of the theory is found to consist 
of the basic solitons, scalar states (or breathers) and excited (or 'breathing') solitons. It 
is also noted that the construction of the S-matrix is valid for any representation of the 
Hecke algebra, allowing the definition of restricted S-matrices, in which case the theory is 
unitary. 
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Introduction 



It is an outstanding problem to understand non-perturbative effects in quantum field 
theory. In two-dimensions the situation promises to be more tractable. Massless Euclidean 
quantum field theories, which describe the critical scaling behaviour of two-dimensional 
statistical systems, exhibit an infinite dimensional symmetry, described by the group of, 
possibly extended, conformal transformations. In a sense, the infinite problem of quantum 
field theory is reduced to a finite problem involving the representation theory of an infinite 
algebra [1] . In the space of massive two-dimensional field theories, there is a class of theories 
that shares the property of integrability with the conformal field theories. In a massive 
integrable theory there exists an infinite number of commuting conserved charges, and 
so there exists some transformation to action-angle variables, and the theory is separable 
(at least classically). For the relativistically invariant classical integrable theories, it is 
interesting to speculate as to the nature of the associated relativistic quantum field theories. 
It is thought that the property of integrability will generally survive quantization and the 
resulting theory will be particularly simple because its S'-matrix will factorize [2,3]. 

For all their simplicity, it has still proved to be very difficult to quantize the classical 
integrable theories. There are a number of approaches, for instance the quantum inverse 
scattering formalism [4] or the formalism of ref. [5] . Another approach, which has proved 
particularly useful, might be called the 'direct S-matrix approach' [3]. The success of this 
latter approach, rests on the fact that the constraints of S'-matrix theory, along with the 
hitherto mentioned property of factorizability, strongly constrain the allowed form of the 
S-matrix, in fact to such an extent that it becomes possible, with a bit of foresight, to 
conjecture a form for the S'-matrix, up to possible ambiguities of CDD type. 

The S'-matrix approach is particularly well adapted to quantizing classical soliton the- 
ories because in such theories there is a direct relationship between the semi-classical limit 
of the S- matrix and the classical soliton scattering solutions [6]. In a classical soliton the- 
ory the scattering of solitons, as described by multi-soliton solutions to classical equations, 
can be described in a very simple way. In the distant past, an iV-soliton solution has the 
form of N well separated solitons of given velocity. In the far future, the solution also 
consists of N solitons with exactly the same distribution of velocities; this is due to the 
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existence of conserved higher spin integrals of motion. If the scattering solution is analysed 
in detail, then one finds that only the centres-of-mass of the solitons change (consistent 
with the overall conservation of momentum). A shift in the centre-of-mass of a soliton 
can alternatively be thought of as a time-delay or advance. The important point is that 
the scattering of the N solitons can be thought of in terms of ^N(N — 1) elementary 
pairwise scatterings, each of which contributes to the shift in the centres-of-mass of the 
participating solitons. 

The picture in the quantum theory is thought to be very similar in the sense that 
the scattering of N solitons can also be thought of in terms of pairwise scattering, with 
the property that individual velocities are conserved (this is the factorizability assumption 
[2,3]). In the quantum theory there can be non-trivial mixings due to mass degeneracies 
(such processes are seen in the classical theory as solutions involving a complex time tra- 
jectory [7]) and the analogue of the time-delay is now played by a non-trivial (momentum 
dependent) phase factor. The semi-classical limit of the phase factor is directly related to 
the classical time-delay. So there is a very concrete connexion between the S'-matrix and 
classical scattering theory. The idea is to use the axioms of S'-matrix theory and the semi- 
classical limit in order to propose a form for the S'-matrix of a classical soliton scattering 
theory. 

The known class of relativistic integrable field theories includes the sine and sinh- 
Gordon theories, the Toda field theories, the chiral models and various fermion models. In 
this work, we will consider the affine Toda theories associated to the algebra s/(n), which 
include the sine and sinh-Gordon theories as special cases when n = 2. 

The sinh and sine-Gordon theories are almost completely understood, in the sense 
that their complete spectra and ^-matrices are known [3,8]. Since we shall argue that the 
more general Toda theories share many features of the sine/sinh-Gordon theory it is worth 
discussing these theories in some detail. The equation of motion for both theories may be 
written 1 

□^ = -^sin(v / 2/30). (1) 

The field (j)(x, t) is a scalar, whilst m and (3 are coupling constants. The sinh-Gordon 
theory differs only in that (3 is taken to be purely imaginary. 

The spectrum of the sinh-Gordon theory is particularly simple, there is only a single 

1 The notation used is not conventional, usually (3 — > /3/y/2, however, it will prove more useful 
for comparing with the sl(n) generalization. 
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scalar particle, which is identified with 4> the fundamental field of the theory. The situation 
for the sine-Gordon theory is, in comparison, much more complicated. Classically the sine- 
Gordon theory admits soliton, or kink, solutions of the form 

0(x,t) = ±H^tan- 1 (e^-*-O), (2) 

where v, a and £ are constants, and cr 2 (l— v 2 ) = Am 2 . The two solutions in (2) represent the 
soliton and anti-soliton. Multiple soliton solutions also exist which describe the scattering 
of solitons. There are other classical solutions, known as the breather or doublet solutions, 
which have the interpretation of a soliton and anti-soliton oscillating about a fixed centre 
with some period. In a remarkable series of papers by various authors [8,9], the complete 
spectrum in the quantum theory has been determined. One has to distinguish two regions; 
the attractive regime which, in our conventions, is < (3 2 < 2ir and the repulsive regime, 
2tv < (3 2 < 4% (the language refers to the remarkable connexion of the sine-Gordon theory 
with the massive Thirring model [10]). The theory is not well-defined when (3 2 > 4tt. In 
the attractive regime, the spectrum consists of a soliton and anti-soliton of mass 

8m 2m 

The first term is the classical soliton mass, whilst the second is the first quantum correction. 
There are arguments to suggest that (3) is actually exact to all order in f3 2 [9]. In addition 
to the solitons there is series of soliton anti-soliton bound states that arise from quantizing 
the breather solutions, yielding a discrete mass spectrum 

m(k) = 2Msin , k = 1, 2, . . . , < ^, (4) 



where 



7 = T^v^- ( 5 ) 



Remarkably, the ground state of the breather system is identified with the fundamental 
particle of the sine-Gordon Lagrangian. Indeed, m(l) = 2m + 0(/3 2 ), where the terms 
0((3 2 ) agree, order by order, with the loop corrections to the mass of the fundamental 
particle. This theory exhibits the phenomenon of nuclear democracy; the solitons can be 
though of as coherent excitations of the fundamental particle, which itself can be though 
of as the bound state of a soliton and anti-soliton. The full 5-matrix for the combined 
breather /soliton system was found in [3]. 
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In the repulsive regime the situation is not so straightforward [8] . The vacuum becomes 
unstable and the spectrum is changed. New scalar states appear which cannot be thought 
of as bound states of the solitons. 

The sine-Gordon theory is also remarkable in that it can be reformulated as the theory 
of a massive interacting fermion, namely the massive Thirring model. In this picture the 
solitons are represented by a Dirac spinor field and the breather states are then fermion 
bound states [10]. 

We now turn our attention to the sl(n) generalizations of the sine/sinh Gordon theo- 
ries. These are the affine sl(n) Toda field theories. The equation of motion of the theory 
is 

2 n 

D^=-^a^> (6) 

P j=i 

The field <f>(x,t) is an n — 1 (= rank s/(n))-dimensional vector. The inner products are 
taken with respect to the Killing form of sl{n) restricted to the Cartan subalgebra. The 
«j's, for j = 1, ... ,n — 1 are the simple roots of sl(n); a n is the extended root (minus 
the highest root) 2 . The fact that the extended root is included in the sum, distinguishes 
the affine theories from the non-affine ones. In the affine theories the ctj's are linearly 
dependent: 

n 

5> i = 0. (7) 

i=i 

Notice that the sinh-Gordon theory is recovered by choosing the algebra to be sl(2). 

For the moment we consider the theories which generalize the sinh-Gordon theory, so 
(3 is real. In the weak coupling limit we expect the theory to include n — 1 scalars with 
masses 

rh a = m a ( 1 + —cot- H ), a = 1, 2, . . . , n - 1, (8) 

\ An n I 

where the first term is the 'classical mass': 

m a = 2msin^— j j = 1, 2, . . . , n — 1. (9) 

An S- matrix for the scalars was conjectured in [11], which is consistent with the Feynman 
diagram expansion in the limit of weak coupling. The S'-matrix only has poles on the 



For convenience we shall think of the labels on the ay's as being defined modulo n, so that 

atj = a j+n 
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'physical strip' corresponding to the exchange of the scalar particles, and so the spectrum 
is complete. This just generalizes the situation for the sinh-Gordon theory. 

Our interest is in soliton theories for which we need to generalize the sine-Gordon 
theory. To do this we can consider the Toda equations when (3 is purely imaginary, so 
(3 = i/3. For n > 2 this means that the equations are now complex. Classically, there 
exists a well-defined set of classical soliton solutions with real masses. Rather than take 
the usual route to quantization we shall proceed directly to an S matrix, showing that 
our proposal is consistent, in the semi-classical limit, with the classical scattering theory. 
Since the Hamiltonian of the theory is a complex quantity (for n > 2), we expect that the 
theory is non-unitary; indeed the issue of unitarity can be addressed directly at the level 
of the S-matrix. 

Our proposal for the S'-matrix involves the intertwiner of the sl(n) quantum group, 
and hence the Hecke algebra. In fact the deformation parameter of the quantum group is 
the coupling constant for the solitons, such that in the weak coupling limit of the solitons 
the quantum group reduces the universal enveloping algebra of the Lie algebra sl(n) and the 
Hecke algebra reduces to the symmetric group, which in physical terms means that the S'- 
matrix is just a permutation and hence describes a non-interacting theory. In crude terms 
the S-matrix has a quantum group symmetry, as well as certain momentum dependent 
symmetries which lead to new integrals of motion. We also point out that the S-matrix 
can be constructed for any representation of the Hecke algebra. 

The plan of the paper is as follows. In §1, we review the construction of the classical 
soliton theory [12]. §2 provides an introduction to the sl(n) quantum group and the Hecke 
algebras, as well as proving a few results which are needed in the following sections. A 
factorizable S-matrix is constructed in §3, which in §4 is proposed as the soliton S-matrix 
of the complex sl(n) Toda theories. Various properties of the S-matrix are investigated, 
for instance we show that the spectrum consists of additional scalar states as well as 
excited solitons. In addition we show, as expected, that the S-matrix describes a non- 
unitary quantum field theory. In §5 we show that the S-matrix is consistent with the 
semi-classical quantization of the classical scattering theory, via a WKB approximation. 
Finally, some comments are made in §6. 
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1. The Classical Solitons 



The soliton solutions of the complex sl(n) Toda theory were constructed in [12]. Here, 
we recall some of the details in order to motivate the construction of the S'-matrix in 
following sections. 

The equation of motion of the complex sl(n) affine Toda field theory is 

2 n 

□^-^E^e'*"-*. (1-1) 

This equation is integrable in the sense that there exists an infinite number of Poisson- 
commuting conserved charges (see for example [13]), one of which is interpreted as the 
energy. This is the integral of the density 

2 n 

H = \ ((d t d>f + (d x ct>f) - ™ £ (e»°r* - 1) . (1.2) 

' 3 = 1 

Notice that the field is periodic with respect to the weight lattice, A*, of sl(n). More 
precisely 

2-7T 

(f>~(f>+—w, Vw G A*. (1.3) 

The constant field configurations <j> = (2ir/0)w, \/w £ A*, have zero energy. There exist 
kink, or soliton, solutions which interpolate between these configurations as x goes from 
— oo and oo. One can associate a topological charge 

t=£- dx d x 4> e A*, (1.4) 

to each solution. 

Explicit expressions for the soliton solutions were found in [12] using the Hirota for- 
malism [14] . The idea is as follows: first one performs a change of variables 

1 n 

% P j=l 

The equation of motion (1.1) is equivalent to 

f jTj ~ r? - t'-Tj + rf = m 2 (tj_iTj +1 - rf), j = 1, . . . , n, (1.6) 
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where, as for the roots, tj = Tj +n . (1.6) is an equation of 'Hirota bilinear type' [14]. The 
characteristic polynomial in this case is 

A, a) = a 2 - A 2 - 4m 2 sin 2 (— ) , (1.7) 

where a and A are continuous variables and oG{l,2,...,n - 1}. It is useful to introduce 
v = A/cr, which is the velocity of a soliton. The iV-soliton solution is then written in terms 
of functions p = 1, 2, . . . , N, one for each soliton, and associated to each soliton 
pair, 1 < p < q < N, where 

( \ 2ttz 
®j{x,t) = a p x - X p t H a p j + i p , j = l,...,n, (1.8) 

subject to the constraint 

T{o v , A p ,a p ) = 0. (1.9) 

In the above £ p is a constant, whose real part represents an arbitrary shift in the centre-of- 
mass, and whose imaginary part will determine the topological charge of the soliton. The 
'interaction function' is given by 

exp7 (p<?) _ -^(op ~ &g, K ~ Aq, a p - a q ) ^ (1.10) 
F{v p + cr q , Xp + X q , a p + a q ) 

where T is the characteristic polynomial in (1.7). Another useful way to write the inter- 
action function is 

sin (0 + *(a -a g ) \ sin (9_ *(a p -a q ) 



/„„\ „ „. \ 2i 2n I \2i 2n 

exp7 (p9) (6) = ) ( ) (- , (1.11) 

where 9 = 6 P — 9 q is the rapidity difference of the two solitons 3 . The general N soliton 
solution is then 

1 1 / n \ 

= E ••• E ex p + E ^ ipq) ■ ( L12 ) 

Mi=0 AtJV=0 \p=l 1<P<Q<N J 

Let us now analyse the one-soliton solution in more detail. The explicit expression for 
the solution is 

0(M) = E a o lo S j 1 + ex P ~vt)+Z+ } , (1-13) 



The rapidity is related to the velocity by v = tanh(9, i.e. tanh^ = A/cr. 



with 



a 2 (l-v 2 ) = Am 2 sin 2 




) 



(1.14) 



Clearly the solution represents a kink, whose centre is at vt — cr _1 Re(£), moving with ve- 
locity v and with characteristic size a -1 . (1.14) simply expresses the relativistic invariance 
of the theory: the the faster a soliton goes the narrower it becomes. 

We now find the topological charge of the elementary soliton, that is the soliton with 
a = 1. Assuming a > 0, as x — > — oo </> — > 0. The limit as a; — > oo is more subtle and 
depends upon the choice for the imaginary part of the constant £. One obtains n different 
values, for instance with the choice 



At this point it is convenient to introduce the weights of the n dimensional representation 
of sl(n), ej (again ej +n = ej). These vectors have the property that Y^=i e j = 0, and the 
simple roots, along with a n , are given by otj = ej — ej+i. By using the n possibilities for £ 
one finds that the topological charges of the elementary soliton fill out the set of weights 
of the the n dimensional representation. 

Similarly one can show that the solitons for other values of a have topological charges 
in the set of weights of the a th fundamental representation of s/(n) 4 . In addition to 
the soliton solutions one can also construct breather solutions. They are obtained by 
considering the two soliton solution, for solitons of equal mass, i.e. either two solitons of 
type a or one of type a and one of type n — a. In the centre-of-mass frame the breather 
solution is obtained by taking Ai = — A 2 = iu), for to e R. The solution is interpreted as 
the two solitons oscillating about a common centre. The new feature of the more general 
theories is that breather solutions can have non-zero topological charge, and so there are 
'breathing soliton' solutions. 

One can easily calculate the masses of the soliton solutions from the Hamiltonian (1.2). 
The resulting masses only depend on the particular fundamental representation that the 
topological charge of the soliton lies in, one finds [12] 




m = 1, 2, . . . , n n odd 
m = 1, 2, . . . , n n even. 



(1.15) 




a = 1, . . . , n — 1. 



(1.16) 



4 The fundamental representations are the irreducible representations whose highest weights 
are dual to the simple roots. 
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Notice that these masses are proportional to the masses of the Toda particles of the real 
coupling constant theory (9) . It seems rather miraculous that the resulting expressions for 
the masses are real considering that the Hamiltonian is in general complex. However, this 
can be traced to the fact that the soliton solutions satisfy a reality condition of the form 
(p* = —M<f>, where M is an orthogonal transformation with M 2 = 1, which acts as a Z 2 
permutation on the roots aj, j = 1, . . . , n. The presence of the involution M implies that 
energy of a soliton is real, since H((j),(3)* = H(—M(j), —(3) = H((p,[3) due to invariance 
of the Hamiltonian under permutation of the roots. Some issues concerning the reality of 
the complex Toda field theories are addressed in ref. [15]. The fact that the masses are 
proportional to the masses of the fundamental Toda particles is explained in ref. [12], via 
a relationship between the bootstrap equations of the scattering theory of the fundamental 
particles and a purely classical analogue of the bootstrap equations for the soliton theory. 

For the case when n = 2 the doublet of soliton solutions of (1.13) reduce to the soliton 
and anti-soliton solutions of the sine-Gordon theory in (2). 

The multi-soliton solutions describe the scattering of the single solitons. For example, 
consider an N soliton solution with v\ > • • • > v^. Initially, i.e. as t — > — oo, the 
solution is approaches a set of isolated solitons, in the order 1, 2, . . . , N as x goes from — oo 
to oo; the position of the p th soliton at time t being 



Finally, i.e. as t — > oo, the solution also approaches that of N isolated solitons moving with 
the same set of velocities as those in the intial state, but now with order N, N — 1, . . . , 1 
as x goes from — oo to oo, the position of the p th soliton being 



If the solution is analysed in detail the following picture emerges. Each soliton retains its 
integrity except when near another soliton, within a distance ~ m . In this interaction 
region the solution is complicated, however, the solitons emerge with the same velocities 
and the only effect of the interaction is to shift the centre-of-mass of each soliton (with 
the overall centre-of-mass being preserved). So an iV soliton scattering solution can be 
analysed in terms of ^N(N — 1) two-soliton scatterings. This is illustrated in figure 1, 




(1.17) 




(1.18) 
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Figure 1. A multi-soliton scattering process. 

where the circles indicate the regions where the solution is not approximated by a set of 
isolated solitons; these are the 'interaction' regions. 

The shift in the centre-of-mass of q th soliton as it 'interacts' with the p th soliton is 
—ry(pq) I Another way to interpret this is to say that the q th soliton on interacting with 
the p th soliton will experience at time-delay given by —^ pq ^/{a q v q ). This shift in the 
center-of-mass is illustrated in figure 2. 

Figure 2. Two soliton scattering 

The question that we now address is: is there a factorizable quantum soliton S-matrix 
whose semi-classical limit yields the time-delays of the classical scattering theory described 
above? In our search for the S-matrix we are guided by the situation for the sine-Gordon 
soliton S-matrix. It is known that the S-matrix in this case is, up to a scalar factor, the 
intertwiner of the s/(2) quantum group. So it appears that we should consider the sl(n) 
quantum groups in order to construct the more general S-matrices. In the next section we 
consider some relevant facts about the sl(n) quantum group and its commutant, the Hecke 
algebra. In §3 and §4 we go on to propose a form for the factorizable S'-matrix describing 
the quantum solitons. 

2. The s/(n)-Quantum Group and Hecke Algebra 

We have seen that classically each soliton mass eigenstate is associated with a funda- 
mental representation of s/(n), in the sense that the allowed topological charge is a weight 
of such a representation. In the quantum theory we expect that the asymptotic state 
representing a soliton to carry two quantum numbers, the velocity (or rapidity) and the 
topological charge. Therefore associated to each external state is a vector in one of the 
fundamental modules of sl(n). We denote the module corresponding to the fundamental 
representation having Dynkin labels 'one' over the a th simple root, and zero elsewhere, as 
V a . The soliton with mass M a is then associated to the module V a . 

Since the theory is integrable we assume that the S'-matrix is factorizable. This means 
that the individual momenta of each external state is separately conserved, so there is no 
particle creation and all the elements of the S-matrix can be deduced from the two-body 
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process [3]. Notice that this is analogous to the classical scattering theory of the solitons 
discussed in the last section. In particular, a quantum scattering process may be analysed 
in terms of pairwise scattering, so the multi-soliton S'-matrix can be constructed in terms 
of the two-body S'-matrix. From a Lie algebraic point of view, the two-body S'-matrix 
must act as the inter twiner: 

S a > b : V a ®V b »V b ®V a . (2.1) 

Notice that if there were no interaction then the ^-matrix would simply permute the vector 
spaces. This gives a clue about how to construct a non-trivial S'-matrix, since there exist 
natural generalizations of the permutation groups known as the Hecke algebras, which 
reduce to the former in some limit. 

In the rest of this section we introduce the Hecke algebras and the associated quantum 
groups. Our approach follows that of M. Jimbo in refs. [16]. 

Central to the subject is the Yang-Baxter Equation (YBE) which can be written as 
an identity in End(V ® V <g) V), for some vector space V: 

(R(x) <g> /)(/ <g> R(xy))(R(y) <g> I) = (I <g> R(y))(R(xy) <g> I) (I <g> R(x)), (2.2) 

where R(x) G End(F <E> V), I is the identity in End(V), and x G C* is the spectral 
parameter. We are interested in the solutions for which R(x) is trignometric in u = \ogx, 
and when V ~ C n is the n dimensional module of the Lie algebra sl(n). The solution 
depends upon an additional parameter g, the deformation parameter, which for the rest 
of this section we take to be generic (not equal to a root of unity) . 

Trignometric solutions of the YBE are naturally described by a quantum group; in this 
case the g-deformation of the universal enveloping algebra of g = sl(n), denoted U q (sl(n)). 
The parameter q sets the degree of deformation and as q — > 1, U q (g) — > U(g). U q (g) is 
endowed with the structure of a Hopf algebra. That is an algebra homomorphism 

A (m) : U q — ► Uf m . (2.3) 

/±{m) d e f] ne g the action of U q (g) on tensor products of the n dimensional representation g : 
U q {g) — >• End(F). For generic q (not equal to a root of unity) the irreducible representations 
of U q (g) are in one-to-one correspondence with those of g. The trignometric solution of 
the YBE commutes with (g x g)(A^ 2 \x)), x G U q (g), and so R(x) lies in the commutant 
of A^(U q (g)). 
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In general, the commutant of A( m ) (U q (g)) is equal to the Hecke algebra H m , which can 
be thought of as a deformation of the symmetric group on m objects, S m . A representation 
of H m is generated by {T a , a = l,...,m - 1}, with T a e End(y® m ), subject to the 
following relations 

(T a -q- 1 )(T a + q)=0 

T a T a+ iT a = T a+1 T a T a+1 (2.4) 

[T a ,T b ]=0 \a-b\>2. 

In the limit q — > 1 we recover the relations for the symmetric group S m , with T a — > a a 
being the generator which permutes the a th and (a + l) th space in the tensor product. In 
general, we can label elements T w e H m with elements uj G <S m , such that T ww > = T W T W > 
(if l{ww') = l(w) + l(w'), where Z(tu) is the length of w £ Sm). 

It can be readily verified that the basic R matrix can be expressed as 

R(x) = xT' 1 - x _1 T, (2.5) 

which satisfies the YBE by virtue of the relations of the Hecke algebra Hz- To make the 
above discussion more explicit we introduce a basis {e^, % = 1, . . . , n} for V. From ref. [16] 
we have 

{q~ 1 e % ®ei i = j 

(<? _1 - q) e i <%> e j + i> j (2.6) 

ej ®ei i < j. 

The representation constructed above actually satisfies an additional constraint, over and 
above the Hecke algebra relations, this is the generalized Temperley-Lieb condition, which 
is the vanishing of the deformed full anti-symmetrizer s~ +1 defined below. 

Just as one can generate higher irreducible representations of sl{n) by considering the 
action of projection operators, formed from elements of the symmetric group 5 m , encoded 
by the Young Tableaux, on the m-fold tensor product of the basic representation, so one 
can form higher irreducible representations of U q (sl(n)) by considering the same process 
with S m replaced by 7i m 5 . In particular, we will be interested in the analogues of the 
fundamental representations g a of U q (sl(n)), a = l,...,n — 1, with g a : U q (sl(n)) — > 
End(V a ) (where gi = g and V\ = V). V a can be projected out of the a-fold tensor product 
V® a with the Hecke algebra analogue of the full anti-symmetrizer: 

V a ~8Z(V* a ). (2.7) 
5 For generic q, not equal to a root of unity. 
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An expression for the Hecke algebra analogue of the full symmetrizer s+ and anti- 
symmetrizer s~ has been given in ref. [16] 

s ± = T Lj2 (±y^ q ^(m-i)/2-i( W )) T ^ (28) 

1 1 weS m 

where [m]! = [m][m — 1] • • • [1] and [m] = (q m — q~ m )/(q — <Z _1 ). are projection 
operators, so that (s^) 2 = s^. For example 

4 = ^(Q + T), s - = ij( ? -i-T), 

and so in terms of these projection operators 

R(x) = (xq - x~ 1 q~ 1 )s^ + (x~ 1 q - xq~ l )s^ ■ (2.9) 

Up till now, we have only considered solutions of the YBE for R(x) E End(Vi ® V\). 
One can also look for solutions to the more general YBE 

(R U2 > Ua (x) ® IurWu* ® R Ul ' U3 (xy))(R u ^(y) ® I Ua ) 

(2.10) 

= (I Us ®R u ^{y)){R u ^{xy)®I U2 ){I Ul ®R u *> Us (x)), 
where R u " u i (x) = a ■ R Ui > u i (x), with 

R Ui ' u '(x) E End(Ui ® L/,). (2.11) 

Here, cr is the permutation a(to (8) u) = u ® w for w E Ui, u E Uj and the Ui are arbitrary 
representations of U q (g). Notice that both the left and right-hand sides of (2.10) map 
Ui ® U 2 ® U 3 to U 3 ®U 2 ® U x . 

Since higher representations can be formed by taking tensor products of the basic 
representation, the higher R w ' ^-matrices can be built out of the basic i?-matrix, via the 
fusion procedure. We briefly describe this process following refs. [16]. 

If we know R u ,u (x) and R u ,u (x) then we can write down a new solution of the 
YBE: 

R u "' u '® u (x) = (I®R u "' u (x yi ))(R u "> u '(xy 2 )®I). (2.12) 

With an appropriate choice of y\ , y 2 we can restrict the above to give a new solution to 
the YBE R u > w (x), where W is in the decomposition of the tensor product U' ® U. The 
choice of y±, y 2 is determined by the requirement that 

W ~ R u '' u (y 2 / yi )(U' ®U)c\U®U', (2.13) 
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i.e. is a proper subspace of U <g) L/ 7 . A simple application of the YBE then shows that 
R u " > w (x) so defined, is indeed a homomorphism from U" ®W to W®U" . So for example, 
notice that since R((— <z) _1 ) oc we have 

V 2 ~R((-q)- 1 )(V 1 ®V 1 ), 

so we can find R}' 2 (x) acting on V\ <S> Vi — (I <S> s^)(V"i <8>Vi<g) Vi): 

R}> 2 {x) = R Vl ' V2 (x) = (I ® R(xyi))(R(xyi(—q)~ 1 ) ® I), (2.14) 
where j/i is arbitrary. In a similar way one finds 

R 2,1 (x) = R V2 ' Vl (x) = (R(x yi ) ®/)(/(g) Rixy^-q)' 1 )), (2.15) 

acting on V2 <S> V\ — (s^ <S> I)(Vi <E> Vi <S> Vi). By repeating the fusion procedure we can find 
R a ' b (x), for a, b = 1, . . . , n — 1, acting in 14 <E> V& ~ (s~ ® ■S5")(^ / i <g> *- a+b - ) ), a subspace of the 
a + b- fold tensor product of V\. 

In general [16], R a,b (x) has the following spectral decomposition 

R a >\x) = R Va > v »(x) = J2pu(x)Pu, (2.16) 

u 

where Pjj is the orthogonal projector, relative to a U q (sl(n)) invariant scalar product, onto 
the irreducible representation U in the tensor product V a <8> Vj, C Vj® ( ^ a+6 - > j and pu(x) is a 
scalar function of a;. Notice that (2.9) is of the form (2.16), because V\ ®V\ = V2© W, since 
s^" is the projector onto the fundamental representation V2, and W ~ s^(Vi ® Vi) is the 
analogue of the symmetric tensor. For generic g, this implies that between fundamental 
representations: 

R a ' b (x) = R v «> v »{x) = p a , b (x)a- +b + ■ ■ ■ , 

where the dots represent non-fundamental representations which appear in the tensor 
product V a <S> Vb, and the full anti-symmetrizer s~ +b is the projector onto the unique 
fundamental representation V a+ b, which appears in the tensor product, (and a + b is taken 
modulo n). 

Notice that at each stage of the fusion procedure we are free to choose an overall shift 
in x, as indicated by the presence of y\ in (2.14) and (2.15). We claim that the following 
application of the fusion procedure defines a solution to the YBE acting in the reducible 
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module V\ © V2 © • • • © V^_i, the sum of the fundamental representations, by recursion 
from the basic solution in equation (2.5): 

R a > b+C (x) = (I <g> R a '\x(-q) c/2 ))(R a ' c (x(-q)- b/2 ) <g> 7) 

(2.17) 

^■"(z) = (7i & ' a (a;(-Q) c/2 ) ®7)(7®7i c ' a (a:(-Q)- 6 / 2 )). 

At each stage we have set the overall multiplicative shift in x. The first equation in (2.17) 
is to be understood as being restricted to V a ®Vb+ c ~ (s~ <g> s b ~ +c )V®^ a+b+C ' > and the second 
to Vb +c <S>V a ~ (^b+c^ s aW^' < ' a+b+C ' > • I n the above we assume that a, 6, c < n and 6+c < n. 
The YBE will be satisfied, by virtue of the fusion procedure, if (2.13) is true, i.e. 



R c ' b 



((-^)- (c+&)/2 )oc, c - +b , (2.18) 



since s b ~ +c is the appropriate projection operator. Since we have not managed to find an 
economical proof of (2.18), the details have been relegated to appendix A. 

In the following section we will need to know the positions of any zeros that R a,b (x) 
might have. Using the fusion equations we have 

#>,i (( _ g) i/ 2) = {R{ _ q) g, J)(J 3 R{1)){s - 7) 
However, 7^(1) = (q — q~ l )I <g) 7 and R(—q) = —(q 2 — Q ,_2 )sJ, therefore 

^,l (( _ ?) l/ 2) = _ {g 2 _ q -2 ){q _ q -l ){s + ^ /)(s - g, j) = o. 

More generally, by using the fusion equations, one can easily show that R a,b (x), for a > 6, 
has zeros at x = (- q )-(a+b-2i-2j)/2 ^ for z = 1, 2, . . . , a - 1 and 7 = 1, 2, . . . , 6. 

Finally we note that using properties of the Hecke algebra alone one can show 

R(x)R(x~ 1 ) = (q~ l x — qx~ 1 )(q~ 1 x~ 1 — qx)I £g> 7, 

and so by using (2.17) we have 

a,b 

R a ' b (x)R b ^ a (x- 1 ) = \ \ g{x{-q)- {a+b - 2l - 23+2)/2 )I®I, (2.19) 

i,J = l 

where 

g{x) = {q~ x x — qx~ 1 )(q~ 1 x~ 1 — qx). 

This concludes our discussion of the quantum group U q (sl(n)) and the Hecke algebra. 
The central result that we will use in the next section is the statement that (2.17) generate 
a consistent solution of the YBE on the reducible module @™Z\V a . 
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3. 



A Quantum Group S-matrix 



An S-matrix constructed from the quantum group intertwiner discussed in the last 
section will form the conjectured S-matrix describing the solitons of the complex Toda field 
theory. The deformation parameter q, in this context, is a coupling constant. The ansatz, 
which was first written down in ref. [17], is a generalization of the soliton S-matrix for the 
sine-Gordon theory, where the appropriate quantum group in that case is the deformation 
of U(sl(2)). For the moment we consider the general problem of constructing an S'-matrix 
out of the sl(n) quantum group intertwiner; the connexion with the Toda theories will be 
addressed in §4. 

In general, the axioms of S'-matrix theory do not strongly constrain the form of the 
S-matrix. The situation in two dimensions for integrable theories is exceptional [3]. To 
start with, integrability implies that there can be no particle creation or destruction, since 
the individual momenta of the incoming particles must be conserved. In addition an 
iV-particle S-matrix can be factorized as a product of ^N(N — 1) 2-particle S-matrices. 
Suppose that the particle states form a number of degenerate multiplets labelled by a set 
of finite dimensional vector spaces {V a , a = 1,2, . . . ,n — 1}, with masses m a . Since the 
masses of the particles must be preserved on scattering the only possible processes involve 
'flavour changing', i.e. the S-matrix is determined by a set of maps 

S a ' b (9) : V a ®Vb — >V h ® V a , (3.1) 

where 9 = 9 a — 9b is the rapidity difference of the incoming particles. Consistency with 
factorization implies that the two-body S-matrix must satisfy the Yang Baxter equation 
(I ® S a ' b (9 1 ))(S a ' c (9 1 + 9 2 ) ® J) (J <g> S b > c (9 2 )) 

(3.2) 

= (S b ' c (9 2 ) ® /)(/ ® S a ' c (9 1 + 9 2 ))(S a > b (9 1 ) I), 
where I is the identity operator on the appropriate vector space. Further conditions come 
from the axioms of S-matrix theory and integrability: 
(i) Unitarity 6 . 

S a ' b (9)S b > a (-9) = I b ®I a , (3.3) 

6 Unitarity in this context does not imply that the underlying quantum field theory itself is 
unitary; that issue depends on the signs of the residues of the S'-matrix at particle poles, an 
issue that we address in S4. 
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where If, <8> I a is the identity in End(Vfe <8> V a ). 

(ii) Crossing symmetry. 

S~ a '\9) = {I® C a ) ■ (a ■ S b > a {m - 9))' 1 ■ a ■ (C- a ® I), (3.4) 

where C a : V a — ► V a , a G {1,2, ...,n — 1}, is the charge conjugation operator, and 
C a C a = I a is the identity operator in V a . As before, u denotes the permutation, 
o~(u <S> v) = v <S> u, and t\ means 'transpose' in the second space which is well defined 
because a ■ S b ' a (9) G End(V b ® V a ). 

(iii) Analyticity. S(9) is a meromorphic function of 9. The only singularities on the physical 
strip, < Im6> < n, are along Re^ = (since there can be no particle production for 
physical values of the rapidity), and the simple poles correspond to direct or cross- 
channel resonances. If S a ' b (9) has a simple pole at 9 = iu c ah in the direct channel we 
say that a particle of mass 

m 2 c = m 2 a + ml + 2m a mb cosu^ b , (3.5) 

is a bound state of a and b. The new particle must itself be included in the particle 
spectrum. If ab — > c can occur then so can ac — > b and be — > a, where bar denotes 
charge conjugation. From (3.5) we deduce the following identity 

<6 + *4 + <c = 2tt. (3.6) 

(iv) The Bootstrap equations. The bootstrap equations give a non-linear relation between 
.S-matrix elements. If S a,b {9) has a direct channel pole at 9 = iu c ah , corresponding to 
a particle in V c = -P^ b (V a ® V&), where is a projection operator, then 

S d < c (#) = (I ® S d ' a (# - mL-))(S d ' b (# + life) ® J), (3.7) 
restricted to <8> V c C <8> <8> V a , and similarly 

5 cd(^ = ®/))(/®5 a ' d (^ + iujg)), (3.8) 

restricted to V c ® Va <Z ® V a ® Va- In the above u c ab = n — u c ah , etc. 

The bootstrap constraints are very powerful because they must be consistent with the 
integrability of the theory. What we mean by this is that the spectrum of possible spins of 
the conserved charges is tightly constrained and this in turn highly constrains the possible 
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masses of the physical states. There is a class of minimal solutions to the above axioms 
for which each particle state is non-degenerate, i.e. V a ~ 1, and the S-matrix has the 
minimum number of poles and zeros needed to satisfy the axioms. Each minimal solution 
is related to a simply-laced Lie algebra. The minimal S'-matrix (when multiplied by a 
function of the coupling constant which introduces no additional poles onto the physical 
strip) is then conjectured to be the S-matrix of the particles of the affine Toda field theories 
for real coupling constant discussed in the introduction. The spins of the conserved charges 
of a theory described by the minimal S-matrix are equal to the exponents of the finite Lie 
algebra g modulo its Coxeter number, and the number of particle states is equal to the 
rank of g. So, for example, the particle spectrum of the sl(n) theory consists of n — 1 
particles with masses 

/7ra\ 

iTia = itio sin y — J , a = 1, . . . , n — 1. (3-9) 

Notice that these masses are, up to an overall scale factor, just the classical masses of the 
fundamental Toda particles for real coupling (9), but also the classical soliton masses in 
the complex Toda theories (1.16). The possible fusions are ab — > (a + b) modn, which 
occur at the rapidity values 9 = iu a b'- 

f— t a + b<n , oinN 

^ = {(2-^)* a + b>n. (3 - 10) 

The charge conjugation operator maps a — > a = n — a. The explicit form for S^ in (9), from 
ref. [11], is 

S a J n (0) = (a + b)(a + b- 2) 2 (a + b - 4) 2 ■ ■ • (|a - b\), (3.11) 
where the following notation has been borrowed from ref. [18]: 



(*) = 



The S'-matrix element S^ in (9) has one direct channel pole at 9 = iu a b corresponding to the 
exchange of the particle a + b mod n, and a cross-channel pole at 9 = iu a i corresponding 
to the exchange of particle a — b modn. Notice that since there is only one possible pole 
in S^ n (9) corresponding to a bound state in the direct channel, we can unambiguously 
write u ab = u c ab . 
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We now construct a series of factorizable S-matrices, with degenerate particle states, 
based on the ansatz 7 

S a > b (9) = X a ' b (9)R a > b (x(9)), (3.12) 

where X a > b {9) is a scalar function of 9, and the R a ' b (x) are the solutions to the YBE 
constructed in the last section. It is convenient to split the scalar prefactor X a ' b (9) into 
two pieces: 

X a, b{e) = S a J Q (9)f a > b (x(9)), (3.13) 

where S^ in (9) is the minimal S'-matrix written down in (3.11). 

The form of the ansatz identifies a particle of mass proportional to m a with the a th 
fundamental representation of U q (sl(n)). We now show that the above ansatz is consistent 
with the axioms of S-matrix theory. One should bear in mind that S^ in (9), by itself, 
satisfies all the axioms of S'-matrix theory. 

The YBE equation for S(9) is satisfied because R(x) itself satisfies the YBE. However, 
we deduce that the spectral parameter x and the rapidity 9 must be related by x = 
exp(c6> + d), for constants c and d. 

Unitarity can be satisfied by an appropriate choice of f a,b (x). Using (2.19) we deduce 
that d = 0, i.e. x = expc6> and 

a, b 

f«*{x) = Yl f(x(-q)- {a+b - 2t - 2j+2 ^ 2 ), (3.14) 

M = l 

where f(x) is a scalar function satisfying 

f(x)f(x- 1 ) = — ^ L— t — (3.15) 

(q 1 x 1 — qx)(q 1 x — qx 1 ) 

The R a,b (x) matrix must be consistent with the bootstrap equations (3.7) and (3.8). 
In addition, the residue of S a,b (9) at the pole 9 = iu a b, corresponding to the direct channel 
bound state c = a + b, must be proportional to the projection operator s~ +b . For example, 
S 1 ' 1 (9) should be proportional to s^" when 9 = iun = 2ni/n. By using equation (2.9), we 
deduce that 

x = (_ g )-»*/2Ti. (3.1 6 ) 

Now we know how x is related to 9 we can rewrite the bootstrap equations in terms of the 
multiplicative variable x. It is easy to show that the two bootstrap equations (3.7) and 



7 We ignore any ambiguities of the CDD type. 
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(3.8) then are identical to the two fusion equations in (2.17). In addition, equation (2.18) 
says that the residue of S a,b (6) at the pole 9 = iu a b is proportional to s~ +b , as required. 

To exhibit crossing symmetry, we have to specify the charge conjugation operator 
C a '■ V a — > V n - a . In fact, V a is naturally dual to V n - a via the action of the Hecke algebra 
analogue of the e-tensor, namely s~. This is because dim(s~(V a <g> V n - a )) = 1, where 
V a <E> V n -a C V® n . From the action T(e^ <g> ej) = ej ® e i? for i < j and equation (2.8) one 
finds 

s~ (ei ® e 2 ® • . • <8> e n ) 



n — n(n— 1)/2 / Q i 7 \ 

9 r 1t V (-(z)'"^®^®...®^ (3 - 17) 

\n ! r ^— ' 



where P n is the set of permutations of {1,2,..., n}, and /(w) is the length of the permu- 
tation w G P n , with respect to simple transpositions. If we define a set of dual vectors 
{e*, i = 1, . . . , n}, with e* ■ = Sij, then the explicit expression for the charge conjugation 
operator is 

C a = Xa E (-9)' ({ij}) ( e i. + i®-"® C 0«® •••»<)• (3- 18 ) 

In the above, A a is a normalization constant determined by C a C a = I a - This gives 

^ jaa^-a(a-l)/2-ii(a-l)/2-aa 



A/J. A/,. — 



\aV.\aV. 



For the case n = 2, corresponding to the sine-Gordon theory, the action of the charge 
conjugation operator is unconventional. In the usual formulation of the soliton S'-matrix in 
sine-Gordon theory, the ansatz of (3.12) needs to be conjugated with momentum dependent 
factors in order to ensure crossing symmetry [19]. However, in the present formulation this 
is not required because the charge conjugation operator acts in an unconventional way: 

Cei = (-q)ie 2 , Ce 2 = (-q)-ie u (3.19) 

to compare with the usual action Ce\ = e 2 and Ce 2 = e\. The unconventional action 
(3.19) makes the introduction of momentum depend factors in the definition of the S- 
matrix unnecessary, but is obviously equivalent to the usual formulation by a redefinition 
of the in and out states. 

To determine the constraint implied by crossing symmetry, it is sufficient to consider 
the relationship between S 1 ' 1 (6) and 5' n ~ 1 ' 1 (6'). It is possible to show directly that the 
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S-matrix elements are related by crossing symmetry if the following equation is satisfied 
by the function f(x): 

n—l n—2 

J] f(x-\-q)^) nO^Hz)"'" 1 - <-ar +1 ) = /(*)• (3.20) 

i=l i=l 

To complete the construction of a consistent S'-matrix, we must find a function f(x), 
which satisfies (3.15) and (3.20). Introducing the notation q = — exp — m\ the explicit 
expression for f(x), derived in appendix B, is 

, M _ r(ffi + A)r(l-*M-A) ~ r(l + isM + (j-l)nA) 

T (ffi + jn\) r (-ffi + ((j - l)n + 1)A) r (l - iffi + (jn - 1)A) 
x r (_^m + jnX ) r (isM + ( (j _ i) n + i) A ) r (l + + (jn - 1)A) ' 

To summarize: the fusion procedure for the .R-matrix provides a solution to the boot- 
strap equations if we identify the particle of mass proportional to m a with the fundamental 
representation g a of U q (sl(n)). 



4. The Soliton S-Matrix 



The quantum group S'-matrix in equation (3.12) represents our ansatz for the soliton- 
soliton S-matrix. We now consider some of the implications of this proposal. First of all, 
the fact that the poles of the minimal part of the S-matrix encode the fusing between 
the solitons implies the quantum soliton masses are proportional to (3.9). This means 
that, up to an overall multiplicative renormalization, the quantum soliton masses M a are 
proportional to the classical soliton masses M a in (1.16), or the masses of the fundamental 
Toda particles m a . 

The topological charge operator is proportional to the Cartan subalgebra generator h 
of U q (sl(n)), which has the following action on V\\ 

hei = eiei, (4.1) 

where is one of the weights of the n dimensional representation. It is immediately 
apparent that 

[h®I + I®h,R(x)] = 0, (4.2) 
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and so the proposed S-matrix conserves topological charge. In the quantum theory there is 
a soliton state for every weight of a fundamental representation. In the classical scattering 
theory, not all the weights of the fundamental representations are obtained (except for the 
n and n representations). 

The fact that the fusing rules of the solitons (3.10) are exactly the same as the those 
of the fundamental Toda particles in the real coupling constant theories, implies that the 
spectrum of conserved charges is the same. All these conserved charges are scalars with 
respect to the quantum group. However, there exist symmetries which are non-trivial 
with respect to the quantum group. These 'residual quantum symmetries' generalize the 
situation for the sine-Gordon theory [19,20]. To see this in more detail we have to introduce 
the generators of U q (sl(n)), {e^, fi, hi, i = 1, . . . , n — l} 8 . (The e^'s should not be confused 
with the weights of the n dimensional representation.) The algebra of the generators may 
be found in ref. [16]. Recall that R(x) is invariant under the action of the quantum group 
whose action is defined by the comultiplication A( 2 ); 

A( 2 )(e,) = q h */ 2 <g> e t + e t <g> q~ h ^ 2 

A (2) {fi)=qhi /2^ fi + f .^ q - hi /2 (4 _ 3) 

A (2) (hi) =hi®I + I®hi. 

So crudely speaking the S-matrix is quantum group invariant. However, the R(x) matrix 
also enjoys a momentum dependent symmetry. To see this we note that the n dimensional 
representation of U q (sl(n)) is identical to that of sl(n). Let (eo, fo) be the raising and 
lowering operators associated to the highest root. R(x) satisfies 

R(x 1 /x 2 )(x 2 1 e q- ho/2 +x 2 2 q ho/2 e ) 

(4.4) 

= (xle <g> q- ho/2 + x\q ho/2 <g> e )R(x 1 /x 2 ), 

where ho = — Ym=i hi- A similar equation holds with eo replaced by fo and x\ and x 2 
interchanged. We interpret the above as a momentum dependent symmetry of the S'- 
matrix. The spin of the generator of such a symmetry follows from the relation between 9 
and x in equation (3.16). With q = exp —iviX, the spin of the generator is nX. 

Since we have the full S'-matrix of solitons we can discover whether there are any 
additional states in the theory, which will manifest themselves as simple poles on the 
physical strip. These new states must then be added into the list of states of the theory. The 

8 Here, hi = a.i ■ h, where h is the Cartan subalgebra generator introduced above. 
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bootstrap equations can then be used to extract the S'-matrix elements of these additional 
states. For instance, in the case of s/(2), the sine-Gordon theory, one finds a set of poles, 
whose position depends on the coupling constant, corresponding to the exchange of scalar 
states. These are the 'breathers' or 'doublets'. The masses of the breathers, as written 
down in (4) , follow from the positions of the poles. So we must search for simple poles on 
the physical strip, over and above the simple poles already interpreted as being due to the 
exchange of solitons. 

To begin with, we consider just the scattering of elementary solitons. There are three 
possible classes of process to consider, (i) The transmission process 9 , ik — > ki, for % ^ k. 
(ii) The identical particle process, ii — > ii. (iii) The reflection process, ik — > ik. These 
processes are illustrated in figure 3. 

Figure 3. Transmission, reflection and identical particle processes. 

The S'-matrix elements for the transmission and identical particle processes, after some 
re-arranging of the arguments of the gamma functions using the product representation of 
the gamma function, are 



00 r (ii. + i _|_ J_) r (ii. _|_ i _|_ t±\ 



X 



V 2tt ~ re ^ nA/ I 2tt n T nA ^ 

r f *i + i + i + x^rf^- + 1 - 1 + i=M 

1 l27r^ 1 ^n^nA>' 1 ^tt^ 1 re ^ n\ ) 



(4.5) 



and 



bii^ii(V) . linX e\ — dik^ki{v)- (4.bJ 

sin 

Consider the S-matrix element Sik^ki(O). It has simple poles on the physical strip (0 < 
lm(0) < tt) at 

e = - 2 -^p+— , p = 0,l,...,[A]. (4.7) 

In the above, the notation [A] means the largest integer less than A. The simple pole at 
9 = 2iri/n, corresponds to the exchange of a soliton with topological charge e$ + ej~ and 
mass M.2 m the direct channel. The residues of the poles in (4.7) are proportional to s% , 
and so they are interpreted as being due to the exchange of excited solitons, or 'breathing 
solitons', in the direct channel with topological charge + e^. The soliton corresponding 



The notation here is short for (ej,efc) — > (efc,ej). 
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to the simple pole at 9 = 2iri/n, is just the ground state of this system. The masses of the 
excited solitons are 

M 2 (p) = 2M X cos (l - I)) , p = 0, 1, . . . , [A] , (4.8) 

where M2 = ^2(0). For the sine-Gordon theory the poles (4.7) and the masses have a 
different interpretation. In this case the states being exchanged are the breathers; their 
topological charge is zero and the masses in (4.8) are equal to the breather masses (4) 
(with the identification 7 = Att/X). 

Similarly the element Sn-,u(9) has simple poles on the physical strip at 

e =-^P> P=l,2,...,[±nA]. (4.9) 

These poles are at x = e~^ ip for which 

R( x = e ~ inp ) = (-) p (q - q~ x )I <g> /, (4.10) 

and so it is natural to interpret these poles as being due to the exchange of scalar states 
{i.e. with zero topological charge) in the cross channel of %% — > ii: these scalar states are 
the analogues of the sine-Gordon breathers. The masses of the states are 

nuip) = 2M X sin (^) , p = 1, 2, . . . , [|nA] , (4.11) 

which generalizes the sine-Gordon result (4). 

We now consider the other S-matrix elements S a ' b (6), where we take a + b < n and 
a > 6, without loss of generality. Rather than give a full discussion, we only present the 
results; a complete analysis will appear elsewhere. There are two types of simple pole, 
those whose positions do not depend on the coupling constant and those whose positions 
do. In the first set are 

7T7 iri 

e = —( a + b), —(a-b), (4.12) 
n n 

which correspond to the soliton a + 6, in the direct channel, and a + 6, in the crossed 
channel. In the second set are 

Z7T7 7T7 

0= -p+-(a + b-2j + 2), p = l,2,...,[±A(a + 6-2j + 2)], (4.13) 

nX n 

with j = 1, 2, . . . , b, and 

Z7T7 7T7 

0=—p+-(a + b-2j), p=l,2,..., [i\(2j-a-b + n)], (4.14) 
nX n 
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also with j = 1, 2, . . . , b. The poles in (4.13) are direct channel poles and those in (4.14) are 
crossed channel poles. The poles correspond to solitons transforming in non-fundamental 
representations, excited solitons in both fundamental and non-fundamental representa- 
tions, and scalar states. The excited solitons corresponding to the fundamental represen- 
tation Qa+b have masses given by the square root of 

M 2 a + Ml + 2M a M b cos (j- (a + b - , p = 1, 2, . . . , [±A(a + b)} , (4.15) 

which generalizes (4.8). The scalar states correspond to simple poles in the direct channel 
of S a,a {9), and have masses 

rh a ( P J)=2M a su>(l(l +j -l)), (4ie) 
j = 1,2,..., a, p= 1,2,..., [l\( n -2j + 2)} , 

which generalizes (4.11). 

Another issue concerns the question as to whether the S'-matrix describes a unitary 
quantum field theory. As we have already pointed out, this is a separate issue from whether 
the S'-matrix is unitary as a matrix. A discussion of this point may be found in ref. [21] 
which discusses the non-unitary field theory describing the Lee- Yang edge singularity. 
Consider the S'-matrix element describing the scattering of two equal mass particles. If S 
has a pole at 9 = iu corresponding to the exchange of particle in the direct channel, then 
for a unitary theory one has 

S W ~ (4.17) 

for some p G R. For a non-unitary theory the residue might have a different sign. For 
example, consider the excited soliton poles of equation (4.7) for the process ik — > hi. The 
explicit expression for the S'-matrix element for this process is equation (4.5). The pole 
corresponding to the p th excited soliton comes from the factor 

iO__l V_ 
2tt n nX 

The only 'dangerous' terms as regards the sign of the residue are 

r| P W-^V.-rf-iY 



n\J \ nX J \ nX J 

which contributes an overall sign of (— ) p to the residue. So the first excited soliton corre- 
sponds to a non-unitary coupling, the second a unitary coupling etc. The message of this 
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result is that the underlying quantum field theory is, in general, non-unitary. This was 
only to be expected, considering the complex form for the Hamiltonian (1.2). However, 
as we have seen for the elementary solitons, there exists a regime for which all the non- 
unitary states decouple. In the next section we make some comments about this for the 
full S'-matrix. 

The S'-matrix reduces to a very simple expression when A = 1. In this case q = 1 and 
so the Hecke algebra reduces to the symmetric group and the S'-matrix becomes trivial in 
the space of topological charges. One can easily verify from the explicit expressions that 



Finally let us consider the case for n = 2, where the S- matrix is the soliton S- matrix 
of the sine-Gordon theory. In this case S m i n (9) has no poles on the physical strip; it is a 
CDD type ambiguity and so for this case alone we drop this part of the S'-matrix without 
affecting any properties of the ansatz. The resulting S'-matrix is exactly that of ref. [3], 
up to the unconventional action of the charge conjugation operator (3.19). 

5. The Semi-Classical Limit 

In this section we will verify that some particular elements of the conjectured S'-matrix 
are consistent with the classical scattering theory. The idea is to use the relation between 
the semi-classical limit of the 5-matrix (the limit as % — > 0) 10 and the time-delays of the 
classical scattering theory. We will only consider the scattering of the elementary solitons 
(the n-dimensional representation) for simplicity. 

As we discussed in §3 there are only three possible classes of process involving elemen- 
tary solitons. The transmission, identical particle and reflection processes. It turns out 
that only the semi-classical limit of the first two processes can be connected in a simple 
way with the classical scattering theory. A discussion of the reflection amplitude requires 
an analysis of a complex time trajectory in the classical theory [7] , which we will postpone 
for a future publication. 

For the transmission and identical particle processes, there is a very simple relation 
between the leading term of the semi-classical limit of the S'-matrix element and the corre- 
sponding time-delay of the classical theory. We simply quote the result which follows from 

10 In this section we shall re-introduce fi into our formulas. 




(4.18) 
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the WKB analysis of ref. [6]. If E is the energy in the channel in question and At(E) is 
the classical time-delay, defining 

/* E 

8(E) = \n B n+\l dE'At(E'), (5.1) 

JE th 

where the number of bound states, or resonances, in the channel is the largest integer less 
than ns (which is denoted [n_e]), and E t h is the threshold energy where the resonances 
are just unbound, then the leading behaviour of the S'-matrix is 

S(6) = exp ^ (8(6) + 0(h)) . (5.2) 

Let us apply (5.1) to the scattering of elementary solitons. In the centre-of-mass the total 
energy is E = 2Micosh(#/2), where 9 is the relative rapidity. The threshold energy is 
2M\. The classical time-delay is from §1 

Writing the integral (5.1) in terms of the rapidity we have 

6(9) = ± nB * -± f de'^ 12 \e'). (5.4) 
P Jo 

In (5.4), ^ l2 \9) is the 'interaction function' (1.11) for the classical soliton scattering. For 
two elementary solitons we have 

^'-^U(^i;^- £) )- 

The two relevant S'-matrix elements are written down in (4.5) and (4.6). In order to 
implement the semi-classical limit we have to know how the coupling constant A (or q) 
depends on fi. We have already seen that the limit q — > 1 can be thought of as the weak 
coupling limit for the solitons, since the S'-matrix becomes trivial, however, this cannot 
correspond to the limit (3 — > because the soliton masses are proportional to f3~ 2 . This 
situation is familiar from the sine-Gordon theory, the point being that the solitons are 
weakly coupled when the coupling constant f3 is large. Mirroring the situation of the 
sine- Gordon theory, we shall find that a behaviour of the form 

\=j^(\o+0(h[3 2 )), (5.6) 
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will be necessary. So as % — > A — > oo. In this limit the product of the gamma functions in 
(4.5) may be approximated by an integral in the exponent, so the leading order behaviour 
is 



Sik^ki(O) -^exp 



n\ 



dx log 



r(f + i + *)r(I + 1 



But this is equal to 



exp 



nAo 
ft/3 2 



r(-j| + £ + *)rH| + 

r(f + l + i + x)r(f + l-l + a;) 

r(f + i + *)r(-f -x)' 



(5.7) 



<ix log 



exp < Z7T- 



ft/3 2 7 
A 



r(g + i-x)r(-^ + x) 

sin {-Tix - y) 
sin (nx - y) 



(5.8) 



+ 



inXr 



de' log 



Kin I t + - J sin ( £ - 



- 2 (|) 



ft/3 2 2yrft/3 2 7 

From, this we can extract the phase shift 

\ _ „ \ f o /'sin (%-. + sin ( - *\ 



2/3 2 ' 47T/3 2 Jo - ~° \ sin 2 (|) 



(5.9) 



Now compare the above with (5.4) and the expression for ^ l2 \9) in (5.5). The two 
expressions agree if 

A = 4tt, (5.10) 

and riB, the parameter which relates to the number of bound states in the channel ik — > ki, 
is equal to A (at this order in /3 2 ). This is consistent with the number of bound states 
found in the direct channel of this process (~ A) found in §4. 

In a similar way we can repeat the analysis for the process ii — > ii. However, it is 
more straightforward to take the semi-classical limit of (4.6) 



= S ik ^ ki (9) exp ^-itt^L^ 



(5.11) 



from which we deduce 



i(0) 



4tt/3 2 



/ de' 

Jo 



log 



\2i n J 



sin 



2 (I) 



(5.12) 
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Again this is the correct semi-classical limit if n^, in this case, is zero. This is consistent 
with the findings of §4, where we found no bound states in the direct channel of the process 
ii — > ii (however, there are bound states in the crossed channel). 

So we have shown, at least for some of the processes, that the conjectured S-matrix 
with the following functional form for q 

q=ew ~w ^ + 0{hp2 ^ ' (5 ' 13) 

does indeed describe the quantization of the classical scattering theory. In order to find the 
higher order terms in (5.13) it would be necessary to go beyond the WKB approximation. 
Such an analysis has been done for the sine-Gordon theory in refs. [9,22] yielding the result 

/ 4n 2 i \ 

q = exp\--^ + Trij , (5.14) 
or equivalently 7 = 4/t/A, where 7 is defined in (5). 

6. Discussion 

We have constructed the classical scattering theory of the solitons of the complex 
affine sl(n) Toda equations. The solitons have topological charges which are weights of the 
fundamental representations of the Lie algebra sl(n), and masses whose ratios are equal 
to those of the conventional Toda particles. Using the sine-Gordon theory as a paradigm, 
we proposed a form for the S'-matrix of the solitons by generalizing the sl(2) quantum 
group to the sl(n) quantum group. The quantum group acts in space of topological charge 
of the quantum soliton theory and the resulting S'-matrix commutes with its action, as 
well as with a momentum dependent operator leading to non-trivial conserved quantities 
over and above the ones associated to a real Toda theory. The ansatz related the coupling 
constant of the Toda theory to the deformation parameter of the quantum group (5.13). 
The quantum masses of the solitons are proportional to their classical masses, up to an 
overall renormalization. This is consistent with the lowest order quantum corrections to the 
soliton masses, which can be computed in the following way. One looks at the linearized 
equation around the soliton solutions. The equation is a multi-component Schrodinger 
equation. Surprisingly, given the fact that the potential is complex, the frequencies of the 
modes are real, and hence the solitons are stable to small perturbations. Furthermore, the 
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zero-point energies of the modes may be summed to give the lowest order correction to the 
soliton masses. Details of this calculation will be presented elsewhere. 

The spectrum of the sl(n) theories, for n > 2, is a good deal more complicated than 
the sine-Gordon theory. In addition to the scalar breather states, there are also excited 
soliton states. Ideally, one would like to construct the S'-matrix elements of these new 
states by applying the bootstrap equations. For the sine-Gordon theory this procedure 
terminates, no new states are then produced, and the complete spectrum just consists of 
the soliton anti-soliton and breathers. The problem of finding the complete spectrum for 
the general soliton theories looks rather formidable. The S'-matrix describes a non-unitary 
theory, as expected from the complex form for the Hamiltonian. For the sine-Gordon 
theory the ground state of the breather is identified with the 'elementary particle' of the 
theory. We now show that such an identification can be made for the sl(n) theories. If we 
expand the equations of motion (1.1) in powers of </> then in the linear approximation we 
expect to see modes corresponding the 'classical' masses of (9), as in the real Toda theories. 
On quantizing we might expect these modes to appear as quanta with mass m a + (9(/3 2 ), 
where the corrections arise from loops. Indeed, we did find scalar particles with masses 
given by (4.16). For each a = 1, 2, . . . , n — 1 there is a discrete spectrum of particles. The 
ground states of these spectra have mass 



In the limit of weak coupling, or H — * 0, using the expression for the classical soliton masses 
and the expression for A we deduce 



which are the masses of the elementary particles in the weak coupling limit. 

From the point of view of the ^-matrix S a ' b (9), there are two types of particle appear- 
ing as bound states, depending on whether the position of the associated pole depends on 
the coupling constant (4.12), or not, (4.13) and (4.14). In the former set, there are only 
the original solitons themselves, associated to the fundamental representations. The latter 
set contains the excited solitons, the solitons corresponding to non-fundamental represen- 
tations and the scalar particles. There is a region for the coupling constant, namely 




(6.1) 




(6.2) 




(6.3) 
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for which all the poles corresponding to the latter states are no longer on the physical 
strip. Remarkably the resulting S'-matrix now describes a unitary quantum field theory, 
since all the non-unitary couplings are associated with the second set of states. Hence, for 
(6.3) the soliton S-matrix is complete and describes a unitary theory. Full details of this 
will be presented elsewhere. 

For the sine-Gordon theory, A < 1 is equivalent to 7 > 4-7T, or f3 2 > 2n, which is 
the repulsive regime, for which the S-matrix, as presented, is no longer valid. It would 
clearly be of interest to discover whether the general Toda theories have an analogue of 
the repulsive regime. 

Returning to the S-matrix constructed in §2, it is easy to see that the axioms of 
S-matrix theory are satisfied algebraically (from the point of view of the Hecke algebra). 
This means that they will hold for any representation of the Hecke algebra. In particular, 
it means that we can formulate the theory in the Interaction Round a Face (IRF), or 
Solid- On- Solid (SOS), picture. From a soliton point of view this corresponds to labelling 
the processes by specifying the vacua between the solitons, rather than the topological 
charges of the solitons. Obviously this is completely equivalent to the 'vertex' point of 
view adopted in §2. However, this equivalence is only true for generic values of q assumed 
in §2. When q is a root of unity, say A = 1/p, where p G Z > n, the 'vertex' description 
is no longer appropriate since some of the elements of the .R-matrix become singular. In 
the 'face' picture, however, there is a way of restricting the allowed variables to a finite set 
for which the .R-matrix is well-defined and the Hecke algebra relations are still satisfied. 
In terms of the solitons this would correspond to restricting the allowed set of vacua to 
some finite set (recall that the allowed set of vacua is isomorphic to the weight lattice of 
sl(n)). For these restricted models the spin of the symmetry generator in (4.4) is equal 
to n/p, hence these theories have fractional symmetries generalizing the situation for the 
restricted sine-Gordon theories in refs. [20]. Notice that the restricted S-matrices would 
lie in the regime (6.3), so they would describe unitary quantum field theories. 

In fact there are more general representations of the Hecke algebra associated to 
certain graphs [23] to which one could also associate a factorizable soliton S'-matrix. In 
addition, one could also consider quantum groups related to other Lie algebras. 
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Appendix A 

In this appendix we prove the result that the R a,b (x) matrices defined in equation 
(2.17) satisfy the YBE. As we have already stated in the text the fusion procedure will 
guarantee that the YBE if (2.18) is satisfied. For reasons of space, the proof that we 
present is not self-contained, since we shall assume two lemmas which may be found in the 
second reference of [16]. 

In the following we use the notation R a (x) = I <g) • • • <g> I <S> R(x) ® / <S> • • ■ <S> I to 
denote the basic R(x) acting between the a th and (a + l) th space in the tensor product 
Vi <S> • ■ • <S> Vi, and s~ b to denote the Hecke algebra analogue of the full anti-symmetrizer 
■s^_ a+1 acting between the a th space and the 6 th space, inclusive, in the tensor product 
Vi <g) • • • <g) Vi. 

Lemma 1. Equation (9) of the second reference of [16]. 
Ra^Ra^ixi-q)- 1 ) ■ ■■R 1 (x(-q)- a+1 )s^ a+1 

=sl a R a (x(-q)- a + 1 )R a . 1 (x(-q)- a+2 )---R l (x), 

and similarly 

R^Riixi-q)- 1 ) ■ ■■R a (x(-q)- a+1 )s^ a 

= S2 - a+lJ R 1 ( a ;(- g )- a + 1 )^ 2 (a;(- ? )- a+2 )... J R a (x). 

Lemma 2. 

= (_g)-a-l _ (_ (? )a+l S 2,a+l^l((-<?) ^^a+l' 

By applying Lemma 2 b times one can easily show 

sr,a + i °c *r, (a.i) 
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and 



We now turn to the main proof. R a,b (x) is defined recursively in terms of R(x) by using 
the fusion equations (2.17). What we need to show is that the R a ' b (x) so defined satisfies 
(2.18): 

R a > b ((-q)- {a+b)/2 ) oc s- a+b . (A.3) 

We proceed by induction. Firstly, equation (2.9) implies that (A.3) is satisfied for a = b = 
1: 

Now consider 

po+1,6 ( (_ Q )-(a+b+l)/2\ 

V ' ) (A.4) 

= (#»' 6 ((- ? )-(°+ 6 >/ 2 ) ® ® J R 1 ' b ((-?)- a - (6+1)/2 ))(^ +1 ® Sb -), 

where we have introduced the projection operator (s~ +1 <g> s b ) to implement the restriction 
V a+ i <S> Vb C ]/ i ®( a+6 + 1 ) j explicitly. Assuming that (A.3) is true for R a ' b (x) we can write 
the right-hand side of (A.4) as 

(s- +b ® J)(J ® J R 1 ' 6 ((-?)- a - (b+1)/2 ))(« a - +1 ® fl6 "). (A.5) 

By applying (2.17) 6—1 times we can express -R 1,& (ir) in terms of the basic R(x) 

R}> b {x) = R b (x(-q)^ 2 )R b ^(x(-q)^ 2 ) ■ • • i^-?)-^ 1 )/ 2 )^. (A.6) 

Using Lemma 1 we can rewrite the right-hand side of (A.6) as 

S - 6 ^ 6 (a;(- 9 )-( b - 1 )/ 2 )^_ 1 (a;(- ? )-( b - 3 )/ 2 ) • • • R^-q)^' 2 ). 

Substituting this into (A.5) we find 

sI, a+b R b U(-q)- a - b )R b+ a-i((-q)- a ^^ 

where we have using the fact that s~ b s~ d = s~ b for a < c < d < b. But by (A.l), the 
corollary of Lemma 2, this is proportional to s~ +b+1 , as required. 
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To complete the proof we must consider 

fta,b+l ( /_ ■ ] -(a+b+l)/2\ 

v H) ) (A.7) 

=(/ ® R^\{-q)-^' 2 )){R a '\{-q)- b -^/ 2 ) ® ® 

The discussion proceeds in the same way. Assuming (A. 3) for R a,b (x) the right-hand side 
of (A.7) is 

(I ® ^ +b )(^ a ' 1 ((-9)- 6 - (a+1)/2 ) ® J)(^ ® s^ +1 ). (A.8) 
By applying (2.17) a — 1 times we have 

£ a >V) = £i(x(-9) ( °- 1)/2 )i22(x(-g) (o - 3)/2 ) • • • £ a (^(-<z)~ (a ~ 1)/2 K,a- 
Using Lemma 1 we can rewrite this as 

Plugging this into (A.8), we have 

^a + 6 + l^((-?)- a - 6 )^((- 9 )- a - b+1 )--^a((- ? )- 6 - 1 ) Sa - +1 , a+6+1 . 

But by (A. 2) this is proportional to s~ +b+1 , as required. This completes the proof. 
Appendix B 



In this appendix we find the function f(x) that satisfies equations (3.15) and (3.20). 
By using (3.15) we may rewrite (3.20) as 

n— 1 n— 1 

i=0 i=l 

Now we introduce u(x) = ((—q)~ 1 x~ 1 — (—q)x)f(x), where (3.15) implies u(x)u(x~ 1 ) = 1. 
u(x) satisfies 

Y\u(x-\-qD = n ^r 1 ;"" 1 ^^ 1 - 

i=0 i=i K HJ v Hl 
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In order to solve this equation for u(x) we introduce the Gamma function representation 
of the sine function. Defining q = — exp — in\ and x = exp— Z7r/x, so fx = —nX6/2m\ we 
have 

x-\- q )- a - X (- q r= 2ni 



r(/i + oA)r(i- fi-axy 

in terms of this 

n"^ 1 , -if ^ TT T(fx + iX)T(l- fx-iX) 

„/ (l { - q) > = nr(-M + (i-OA)r(i + M-(i-<W 

To solve this equation we introduce the notation 

n-1 



i=0 

in terms of which 

11 vp 1 = 1 wQx + jnX) 

and similarly for the other terms in (B.l). Using these facts we can write down 
an expression for u(x) as an infinite product of Gamma functions. In fact because 
((—q)~ 1 x~ 1 — (—q)x)~ 1 = r(/x + A)r(l — fx — X)/2ni the explicit expression for f(x) 
is 

r(^ + A)r(i - fx - x) r(i + » + (j - i)nX) 

J[X> 2m *-\r(l-fx+(j-l)nX) 

3 =i v p u ; ; (B.2) 

( £Qx + jnX)T(-fx + ((j - l)n + l)A)r(l -jx+ (jn - 1)A) 

r(-/x + jraA)r(/z + ((j - l)n + l)A)r(l + /x + (jn - 1)A) - 
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